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Abstract 

We study the quark confinement problem in 2+1 dimensional pure Yang-Mills 
theory using euclidean instanton methods. The instantons are regularized and 
dressed Wu-Yang monopoles. The dressing of a monopole is due to the mean 
field of the rest of the monopoles. We argue that such configurations are stable 
to small perturbations unlike the case of singular, undressed monopoles. Using 
exact non-perturbative results for the 3-dim. Coulomb gas, where Debye 
screening holds for arbitrarily low temperatures, we show in a self-consistent 
way that a mass gap is dynamically generated in the gauge theory. The mass 
gap also determines the size of the monopoles. In a sense the pure Yang- 
Mills theory generates a dynamical Higgs effect. We also identify the disorder 

operator of the model in terms of the Sine-Gordon field of the Coulomb gas. 
PACS numbers: 12.38 Aw, 11.15 Tk 
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I. INTRODUCTION 



The problem of quark confinement is one of the "old" unsolved problems in theoreti- 
cal physics. Despite intense activity over the past two decades, and several approaches to 
the problem, it is surprising how little we know about this phenomenon. Lattice gauge 
theories, together with the theory of renormalization group (which provides the basic con- 
ceptual framework for all other approaches as well) is the only known quantitative and 
reliable method of attacking this problem and Monte Carlo simulations have indeed almost 
"demonstrated" confinement in pure gauge theories in four dimensions. However there is the 
nagging feeling that we can "demonstrate" confinement, even calculate relevant quantities 
to some degree of accuracy - but we still dont "understand" confinement. Despite impressive 
progress lattice methods remain a black box. More specifically we do not understand in full 
generality whether (if any) some specific type of gauge field configurations are responsible 
for confinement and whether we can arrive at a consistent picture of the vacuum of strongly 
coupled gauge theories. 

Physical mechanisms for confinement are, of course, as old as the idea itself. One of 
the most significant ideas, proposed by Mandelstam 't Hooft and Nambu || is that 
of dual superconductivity - that the gauge theory vacuum is a condensate of magnetic 
monopoles. By the dual version of Meissner effect, quarks would then be naturally confined. 
The formulation of this idea gave rise to deep results concerning duality between electric 
and magnetic fluxes and its implications to the phase diagram of general gauge theories [f§] . 
The idea is very convincing, and lattice results indeed seem to support it : but is rather 
difficult to establish in usual gauge theories. Recently, however, there has been significant 
progress in supersymmetric gauge theories where holomorphy and global properties of the 
moduli space of vacua have been used to argue for non-perturbative phenomena like quark 
confinement and chiral symmetry breaking ||. This approach is essentially Hamiltonian : 
one tries to obtain a picture of the vacuum wave functional. 

A complementary viewpoint stems from the Euclidean approach to the problem. This 
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is the idea that Euclidean instantons essentially disorder the vacuum and lead to color 
confinement. The only known successful implementation of this idea is the classic work 
of Polyakov |J who showed that confinement indeed occurs by this mechanism in three 
dimensional SU(2) Yang-Mills theory coupled to adjoint representation Higgs field. The 
Higgs breaks the gauge group to U(l). The instantons of this model are the 't Hooft- 
Polyakov monopoles. (In the three dimensional Euclidean theory these are of course not 
solitons, but tunnelling configurations). When the mass scale of this symmetry breaking - 
the mass of the W boson, mw is large compared to the mass scale set by the dimensional 
gauge coupling, a dilute gas of monopoles provides a self-consistent picture. The resulting 
monopole plasma leads to Debye screening. Wilson loops in the fundamental representation 
obey an area law and a careful treatment shows that the adjoint representation Wilson loops 
obey a perimeter law [[?]]. This is exactly what one expects. The argument can be extended 
to the SU(N) theory as well. 

In this paper we extend this approach to pure Yang-Mills theory. This is a much more 
nontrivial system for several reasons. Firstly, perturbation theory, though ultraviolet finite, 
is hopelessly infrared divergent. Secondly, because this theory has only one length scale set 
by the gauge coupling g, one does not have the luxury of having another length scale m-w 
to enable a controlled semiclassical approximation. The classical monopole configurations 
are singular configurations in the continuum limit and hence the renormalization of the 
monopole gas is very nontrivial. 

It is possible to regulate the singularity by modifying the fields inside a "core" of some 
size A. The classical action now depends on A. We do not vary the size A of the monopoles. 
Rather we treat it as a parameter of the theory to be self consistently determined by the 
mass gap. 

If the fluctuations around a regularised monopole solution are decomposed in terms of 
representations of the direct product of spin and color groups, the even parity S-wave fluc- 
tuations are unstable ||. This is because the background magnetic field has a long range 
Coulomb tail. However the fluctuation problem in the Yang-Mills theory should be per- 
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formed not around a single monopole solution, but around the neutral plasma of monopoles 
which populate the vacuum. These monopoles have long range Coulomb interactions : other 
monopoles affect the field around a typical monopole in a nontrivial way. More significantly 
the monopole positions and charges are fluctuating, which make the charge density field of 
the monopoles a dynamical variable. Such a fluctuation problem is difficult to solve exactly. 

In this work we address this question in an approximation guided by the physics of the 
problem. We incorporate the fluctuations of the charge density p(x) by invoking known 
rigorous results of the three dimensional Coulomb gas due to Brydges §| . In [|] it has been 
shown that for a given arbitrarily low temperature there is a chemical potential (fugacity) 
such that the correlation of the density operator cluster. In other words there is Debye 
screening at all temperatures. The field around a given charge thus decays exponentially 
over a debye length l D as opposed to a power law decay. 

Our strategy is as follows. The aim is to show that Debye screening is self consistently 
realized in the plasma of magnetic monopoles. The main complication is that the fugacity 
of the Coulomb gas is itself a functional of the density of monopoles p(x). However using 
the results of || we can argue that for a given value g 2 of the gauge coupling there exists 
a fugacity of the monopoles for which the plasma has a finite Debye length. The mean 
'magnetic field' is also screened with a fall off given by the Debye length. In this sense a 
monople configuration in a plasma does not have a Coulomb tail and such a configuration 
which incorporates this collective property we call a 'dressed' monopole. Small fluctuations 
of the gauge field around a dressed monopole are expected to be stable for reasons similar 
to the Yang-Mills-Higgs theory. There the potential appearing in the stability operator 
decays exponentially with a scale ^— due to contribution from the Higgs field. Hence in the 
pure Yang-Mills theory one seems to have a dynamical Higgs effect that is produced by the 
monopole plasma. We also indicate how the mass gap determines the size of the monopoles. 

Finally we discuss the relation of our approach with the work of 't Hooft |1| . We give an 
explicit representation of the disorder operators in 2 + 1 dimensional Yang-Mills theory and 
indicate that the dual theory is a Z 2 non-linear sigma model. 
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Our discussion of the confinement problem gives a picture of the dominant configurations 
in the Euclidean framework. Feynman [T(| has given qualitative arguments for the ground 
state wave function of this gauge theory in analogy with his work on the roton spectrum in 
liquid helium. It would be interesting to relate these two approaches. 

The plan of this paper is as follows. In Section II we discuss the topology of the vacuum in 
both the Georgi-Glashow model and the pure gauge theory. In Section III we briefly review 
the generation of mass gap and confinement of A-ality in the Yang-Mills Higgs system. In 
Section IV we discuss self-consistent Debye screening of the monopole plasma in the pure 
Yang-Mills theory. Section V is devoted to the dual representation in terms of disorder 
operators. Section VI is devoted to conclusions. In Appendix I we state the main result of 

ti- 
ll. TOPOLOGY OF THE GAUGE CONDITION 

In this section we shall review the topological properties of the vacuum of 577(2) Yang- 
Mills theory in 2 + 1 dimensions following flIE| . The physical degrees of freedom are most 
transparent in a unitary type of gauge. For pure gauge theories this is defined as follows. 
Consider some local operator X(x) which transforms according to the adjoint representation 
of SU(2). The unitary gauge is now defined by 

[X(x),r 3 ] = (1) 

(t 1 ,i = 1,2,3 stands for the Pauli matrices). This gauge condition retains the £7(1) gener- 
ated by r 3 as an unbroken symmetry. If the model contains adjoint Higgs fields in addition 
to the gluons, X may be the Higgs field itself and (p]) is the conventional unitary gauge in 
such a Georgi-Glashow model. In pure Yang Mills theory X has to be constructed out of 
the gauge fields alone. 

We may write the matrix operator X in the form 

X = XI + j2e a (x)r a (2) 

a=l 

5 



where I is the identity matrix. Then the points xq where 



e a (x ) = 



(3) 



are singularities of the gauge condition ([!]). As shown in |13| these singularities are nothing 



but magnetic monopoles with respect to the unbroken U(l). In 3 + 1 dimensions these 
monopoles are "particles" which respresent dynamical degrees of freedom (in this gauge) 
other than the conventional fields. In 2 + 1 dimensions these monopoles are instanton 
configurations which will play a crucial role in determining the properties of the vacuum. 

In the Georgi-Glashow model the monopoles are given by the well known 't Hooft- 
Polyakov solution, while for pure gauge theory they are described by the Wu-Yang solution. 

III. CONFINEMENT IN 2+1 DIMENSIONAL GEORGI GLASHOW MODEL 

We now briefly discuss the salient features of the mechanism of quark confinement in the 
Georgi-Glashow model (for details see Ref. || and 0). 
The model is described by the lagrangian density 



where denotes the gauge field and $ a Higgs field in the adjoint representation of SU(2). 
The instantons of this model are 't Hooft-Polyakov monopoles. For reasons to be discussed 
later we are interested in monopoles with minimal charge :q = ±1 These are the 't Hooft- 
Polyakov monopoles. In the unitary gauge [<&, T3] = the fields are given by: 
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where q = ±1 and 
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—6 sin 6 + 9 cos * 
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- tan - = L> M 

7 Zi 



(6) 
(7) 



where (r, 9, 0) denote polar coordinates in space-times. The functions K(£), H(£) obey well 
known differential equations for (|) - (|6[) to be a classical Euclidean solution. For r 3> m^ 1 
one has K(rm\y) — and H{vr) ~ w. m^/ denotes the "size" of the monopole. 

Note that given a configuration of monopoles and anti-monopoles, the Weyl group 
changes each monopole into an anti-monopole and vice-versa. In principle one may fix 
the gauge further to remove this discrete degeneracy. We, however, prefer not to do so and 
average over the Weyl group. Then one may freely perform a sum over all the q — ±1. 

As shown in , in the dilute gas approximation the path integral may written as a grand 
canonical partition function of a gas of monopoles: 

OO 1 



N=0 JV • 



N 



(8) 



where 



Qn = J2 I II dx * ex P f ~~r 12 1 9a9b 



and the fugacity J is: 



J 
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detD 2 



-1/2 



X 6 



det A FP " 



(9) 



(10) 



^7T " ydet-<9 2 / \det-9 2 / 

D 2 denotes the small fluctuation operator around single monopole field, App is the Fadeev- 
Popov operator, and s is the action of a single monopole. 

In (|^) - flrop we have treated the various zero modes of D 2 by the standard procedure of 
collective coordinates [13]. The Weyl degeneracy of the unitary gauge condition has been 
accounted for by summing over q a = 0, ±1 for each space time point, as noted earlier. 

The Coulomb gas of eqns. (|8|-|T0D may be expressed as a massive scalar field theory using 
the sine-Gordon transform: 

,2 



T>x(x) exp 
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32n- 



J rf 3 x{(V X ) 2 -2M 2 (l-cos X )} 



(11) 



where 



M 2 = — . 12 

T 

Let us now examine self-consistency of the above scheme. When J is small, so that the 



scalar field theory is weakly coupled, M is precisely the mass gap. Equation ( |TT| ) expresses 
the long distance behavior of the theory in a dual representation. The classical action of a 
single monopole has the form 



9 2 



(13) 



When mw/g 2 ^ 1, it follows from (10) that J/g 6 <C 1. Further the "classical" piece in 



J dominates over the contribution from one-loop fluctuations. Thus the scalar field theory 
(pT|) is indeed weakly coupled; M is the mass gap and from (0) M/g 2 is small. 

Since M is the mass gap, the correlation length £ = 1/M. Also, J denotes the probability 
of occurence of a single monopole and hence the number of monopoles in a Debye volume 
of size £, Ng is given by 

»t = W~Ji >l (14) 

where we have used fli"2]). 

We immediately see that there is a large number of monopoles in one Debye volume 
and hence the potential x may be treated classically. The situation here is identical to 
the Debye-Huckel theory of electrolytes in the limit of high temperatures where the Debye 
length is large and the smoothly varying potential field \ satisfies the classical sine-Gordon 
equation. 

The Wilson loop average is given by 

(W(c)) = ^(TrPexpfifAndx^) (15) 

where = A a ^- when the external quarks are in the fundamental representation of SU{2) 
and A^ = A^T a when the quarks are in the adjoint representation. Here T a are the genera- 
tors of the adjoint representation. 
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The "classical" contribution from the Wilson loop factors out: 

(W(c)) = (W{c)) c ^{c)} qu (16) 

where {W{c)) cl = \{Tr Pexp (i f Affidx^), the average being evaluated in the ensemble 
given by (|TT|) . [c<j(c)]g U obeys a perimeter law: 

[u{c)] qu ~ exp(-aP) (17) 

where a is a constant and P = perimeter of loop. 

Let us first discuss the Wilson loop in the fundamental representation. In the ensemble 
given by ( |TT| ) one has 

(W(c))T = J V x (y) expf-3^: / A(V(x - |) 2 - 2(1 - cos*))] (18) 

Where rj(y) is the magnetic scalar potential due to a dipole layer of unit strength on sheet 
S, which is the solid angle subtended by the loop at the point y In (|T8"D we have scaled the 
distances: y = Mx. 

Since M is small, ([18]) may be evaluated by stationary phase approximation. The result 

is 

(W(c))^ d = exp(-aA) (19) 



where 



g 2 M 

a 



327T 2 J- 

and A = area of the loop in physical units.. Then, using (|T6| ) it is easily seen that the 
Wilson loop obeys an area law, indicating that quarks in the fundamental representation 
are confined. 

As shown in |7j the Wilson loop in the adjoint representation obeys a perimeter law 
instead. This happens because there is no distinction between the quantum of magnetic flux 
created by the Wilson loop and that created by the monopoles in the vacuum. The picture 
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of confinement discussed in this section is valid for any SU(N) gauge group. The mechanism 
is essentially Debye screening in a gas of "non-Abelian" monopoles belonging to the adjoint 
representation of a dual group *SU(N) J/J. 



IV. CONFINEMENT IN PURE YANG-MILLS THEORY 

In pure Yang Mills theory there is no Higgs field. This has two consequences. First, in 
the absence of the second scale (i.e. mw), the perturbation expansion is infrared divergent 
and the theory cannot be defined perturbatively in the infinite volume limit. Secondly the 
classical monopole solutions are Wu-Yang monopoles |T5[ which have zero size and infinite 



action. We will regulate these monopoles by assigning a size A, which is explained in the next 
subsection. We will then construct an expansion around a plasma of such monopoles and, 
as explained in the introduction, argue that these monopoles are stable against fluctuations. 



A. The monopole solution 



We will consider monopoles which have some size A. The field due to single monopole 
at x = is given in the unitary gauge by: 



i 

2 



K(r/X) 



r 

K(r/X) 
r 

i e 

— tan - 

r 2 



• cos + 9 sin i 
» sin (f) + 9 cos ( 
>L = D„ 



(21) 



in spherical coordinates. Here K(r/X) is a structure function regulating the fields at r = 0. 
The function K(r/X) goes to 1 as r — > as follows 



Kir IX) ~ 1 - — for r -> (22) 
while at r = A, K(r/X) = K'(r/X) = and remains zero for r > A. Furthermore K'(r/X) is 



continuous at r = A. As shown by Banks, Myerson and Kogut [Tj|, one can choose such a 
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K(r/X) so that the configuration (21) is a classical solution. The action of single monopole 
is 

2 



An f°° , 
s — — I dr 



g 2 Jo 



dK\ / if 2 — 1 
dr J \ 2r 2 



(23) 



Note that the monopole field is abelian outside the monopole core. Consider a gas of 
such monopoles such that if the positions of the monopoles etc., one always has 

|x a — > 2A. The field configuration in the regions outside the core of the monopoles is 
approximated by 

N 



<tf(x) = £ A^(x - x a ) 



(24) 



a=l 



where x a denotes a monopole position and iV is the number of monopoles. In our self- 
consistent approach, we assume that (|24j ) represents the dominant field configuration in the 
euclidean path integral. The total action of this gas is 



S, 



cl 



QaQb 



9 2 a^b 



Finally we record the single monopole field configuration in the radial gauge 



(25) 



Ad 



X J 



[1-K(r/X)) 



(26) 



B. The path integral 



The monopole configuration is used to evaluate the path integral by the saddle point 
method. One expands the field around A C J; which is a classical solution outside the core. 
The form of the solution inside the core is unimportant for our purposes: 



(27) 



The path integral may be formally written as: 



Z = J Y[dAfi(x) exp(—Sc£) exp ^— J aD 2 ad 3 x 



(28) 
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where S c e is given by (p5| ). Here D 2 is the stability operator for small fluctuations: 

J aD 2 ad 3 x = J C"{A cl ) a 2 d 3 x (29) 

where C"(A c e) is the second functional derivative of the Lagrangian density evaluated at 
A c e- 

Equation (28) as it stands is meaningless since D 2 has a zero eigenvalue for each symmetry 
of the original Lagrangian. The local gauge symmetry is fixed by requiring the fluctuations 
to satisfy the background gauge condition 

V M (.A rf K(a;) = (30) 

where V /J (^4. cf ) is the covariant derivative evaluated at the configuration A . 

V,(Aa) = d, + i[Af, ] (31) 

This gives rise to the usual Fadeev- Popov determinant (det App). The zero modes arising 
from the breaking of the global translation invariance and the U(l) transformations obeying 
the background gauge condition but nonvanishing at infinity are replaced by integration over 
corresponding collective coordinates. Finally we have to sum over all N with the standard 
division by N\. One finally has the formal expressions 

OO 1 

Qn=J2 (-7= s 3/2 ) [ n dx a e- s °* J (32) 

!aa\ VV^ J J a=l 



det£> 2 \ 1/2 / det A^p 



J ~ 1 det -d 2 det -d 2 (33) 



C. Instability of the Single Undressed Monopole 



In the usual semiclassical method the dilute instanton gas [I7| is noninteracting and one 
writes 



12 



det D 2 = (det dT (34) 

where d 2 is the stability operator for the single monopole configuration. 

We will see, however, that the d 2 has negative eigenvalues signifying the instability of a 
single monopole. In the background gauge (0) the operator d 2 has the form: 

ci 



d 2 = ^V Q (Ace) V° (Aaj+i F±, (35) 



where A C £ is defined in ([H]) and Fff v is the corresponding field. The Wu-Yang case corre- 



sponds to K = in fl2ip . Following Yoneya || we shall use a spherical basis in the product 
space [(space-time) ® (isospin space)]. In this product space A^ and the fluctuations a M 
become tensors. Then the unstable modes of d 2 are given by 

(a) Odd parity S waves: 



t>i = a ri 



V2 



+ 



(b) Even parity S wave: 



1 

V2 



The tensor indices on refer to the abovementioned spherical basis. The corresponding 
eigenvalue equations are: (See Yoneya, Ref. ||) 

( d 2 3K 2 - 1\ / - n 2/ - . 

+ 5— (r04) = " e 2 (r0 4 ) (36) 



dr 2 r 



2K 2 l r (K d K-tK'\ l 2j 

— r fa-V2[ — - r + 2 04 = ag0! (37) 

r z \ r dr r z ) 



r-(\ d , r , x K-rK'\ , / d 2 2 d if 2 - l\ , 9 , 

v/2 - - rlf) ^ + _ + 4 = a 4 (38) 

\r z dr r J \ dr z r dr r z I 
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^ + -\ <P, = V2-<p A (39) 
dr r r 



For the Wu-Yang monopole K — and the instability is obvious from flHop and (|3"8|). In this 
case both rcf)^ and r04 obey the equation 

( d2 1 \ 2 

1 ~ ^2 J V = <*W, lf) = J-04, ^04- (40) 

The above equation is the Schrodinger equation for a particle in a spherical potential Uo(r) = 
— ^2 with 7 = 1. For such a potential it is known that when 7 > 1/4 there are bound states 
ill 

When the structure function K{r/X) 7^ the odd parity S'-waves become stable. For- 



mally one can invert (|39f) to express 04 in terms of 0i and insert this into ( p7|) to obtain: 



dr 2 r dr j 



In the corresponding Schrodinger problem the "potential" is now positive everywhere and 
hence a 2 > 0. Since one can now rewrite ( |3"9"D as 

, V2 



[ r yK(y)My)dy (42) 

JO 



r z Jo 

stability of 04 also follows. 

However, regularising the field near r = does not remove the even parity S-wave 
instability. In fact (J5BJ) has a potential that becomes asymptotically — -y as r — > 00. For 
such a potential there are still an infinite number of negative eigenvalues. This result is 
independent of the detailed form of the potential near r = 0. Hence the single regularised 
monopole is unstable and cannot be treated as a dominant configuration in the path integral. 

It should be noted in the presence of Higgs fields the potential is replaced by 

3K 2 + H 2 -1 



(43) 



where H(vr) is the function in (j^) . As r — > 00, H{yr) — vr + e~ mHT (where uih = -^^^ 
is the Higgs mass) , and hence the X Coulomb tail is cancelled and we have a screened 
potential e " , which removes the potential instability. 

14 



D. Debye Screening in the Monopole Gas 



Our main point is that when the instantons are interacting, the fluctuation problem can- 
not in general be factored into N copies of the fluctuation problem for an isolated monopole. 
Rather one should consider the stability of the neutral plasma of monopoles as a whole. This 
statement also applies to the Yang-Mills-Higgs system in the previous section. The main 
reason behind this is that we have an integration over the positions of the monopole and 
the charges, or equivalently a functional integration over the charge density field. The effect 
of this averaging over the charge density is very nontrivial. The results of [[| show that 
the charge density field clusters so that the theory of the density field generates a mass gap 
dynamically. We summarize these results in Appendix I. 

In fact the results of || mean that in the neutral plasma the fluctuations of the magnetic 
field are bounded and the \ magnetic field of a single monopole is Debye screened to 6 ™ r , 
where M(g 2 , A, z) is the non-perturbative mass gap, which depends on the coupling g 2 , the 
monopole size A which is effectively the cut-off of the Coulomb gas and the fugacity z. 
Recall that the source of instability for a single isolated monopole is the long range tail of 
the Coulomb potential. One might, therefore expect that in the screened neutral plasma a 
"dressed" monopole whose Coulomb tail has been screened can in fact be stable. 

In the following we shall assume that the monopole gas is dilute. Using translation 
invariance we focus on one monopole at x = x a and its neighbourhood. We are thus 
considering the problem in the presence of a single source at x = x a . Recall that the fields 
outside the monopole cores of size A are abelian. Thus, in the unitary gauge it follows from 
(p4j) and fl2lD that outside the core of this monopole the field is 



where we have set xn = x a . Inside the core the effect of the core field of the other monopoles 
can be ignored and we have 




(44) 
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i / qD^x - x a ) W (x - x a ) 

+ - _ ) (45) 

2 \W+(x-x a ) -qDJx-x a 



where = A 1 ± iqA 2 and A 1 , A 2 are as in (21). Introducing the charge density 



A? 



p{x) = ^q a 5(x - x a ) (46) 



a=l 



and assuming that in the mean, for large N, Pn — Pn-i we can rewrite (fj5|) as 

r ■ 

A^x; x a , \p})= j d 3 y — p(y)D^(x - y) 

( D^X-Xa) W~(x-X a ) \ 

+ 6{\-\x-x a \)\ (47) 

\W+{X-X a ) -D^X-Xa)/ 

In the above expression the sharp 9 function may be replaced by a smoother version. 

Debye screening means that in the presence of a source the density p(y) has a mean value 
p(y,x a ), fluctuations around which are small. A^(x;x a ,[p]) then represents a "dressed" 
monopole configuration. In our case this "source" is provided by the particular monopole at 
x = x a in the plasma and the statement pertains to the field in the neighbourhood of this 
particular monopole. 

The crucial point is that since Debye screening holds, we can assume self- consistently 
that the gas of "dressed" monopoles is weakly interacting, unlike the "bare" monopoles. The 
field around a dressed monopole decay exponentially over a distance scale set by the Debye 
screening length In = -g. If the average distance between the monopoles is much larger than 
Id then the interaction between such dressed monopoles vanishes and the operator D 2 has a 
N-fold degeneracy. In other words the potential appearing in the stability equation resembles 
N far separated potential wells. In this situation we have, using translation invariance, 

detD 2 ~ (detD 2 [p}) N (48) 

where D 2 [p] now denotes the stability operator for a single dressed monopole which is the 
same for any monopole in the plasma. For finite distances between monopoles the exact 
degeneracy is lifted and eigenvalues of D 2 organize themselves in bands. This would lead to 
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corrections to the result (HI) which may be expanded in powers of f 2 - where l m denotes the 

1 1 'm 

average distance between the monopoles. 

In this regard there is a difference between the pure Yang Mills system and the Yang- 
Mills-Higgs system in the Higgs phase. As our stability analysis in the previous subsection 
indicates (see e.g. equation (|43"D) the presence of the Higgs field means that the potential 
appearing in the stability operator around a single monopole decays over length scales of the 
order of Since the debye length is much larger than corrections to the extensivity 
of the small fluctuation determinant appear as powers of ; m . 

E. The Sine-Gordon Transform and Dynamical generation of Mass Gap 

We can now rewrite the theory in terms of a sine-Gordon model. The path integral is 
written as 



e-f Ell/ **«P H^lT^h) ( Q W N ( 49 ) 



where we have defined 0[p] (functional of the charge densities) as 

_ ( det D"[p} y 1/2 ( detA FP (p) \ ^ 

and s is the action for a single monopole. 

It may be noted that we are dealing with a superrenormalizable theory. Thus the ex- 
pression ( |50|) is ultraviolet finite. 

We then have 

N=0 JV ' { qa }i=l J \ it l^a ^b\ / 

where the mean fugacity J is given by the formula 

J = ^= e~~ s ®[p\ (52) 

'7T 



We can now perform the sine- Gordon transform as in section 3, and we write fl5lD as 
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T>x(x) exp 



327T 1 



J d 3 x{(V X ) 2 -2M 2 (l-cos X )} 



(53) 



where 



M 



16vr 2 J 



(54) 



It is important to emphasize in accordance with the discussion of || , that the quadratic 
term / d 3 x(Vx) 2 in (H) must be understood in a regularised sense, so that the Coulomb 
potential between the monopoles is valid only upto distances greater than the core size A. 
In this sense A is the cut off (lattice spacing) of the sine-Gordon theory. 



F. Stability of the dressed monopole 

We now discuss, in some more detail than previously, whether the function Q[p] which 
is used in the definition of the fugacity J in fl52|) is well defined. This issue is important 
because we have already indicated in Section 4 that as we average over p(x) the operator 
D 2 [p] has negative eigenvalues when p corresponds to a single isolated monopole. 

We will now argue that D 2 [p] is a positive operator. Recall the form of D 2 [p] in the 
unitary gauge, outside the core of the dressed monopole which we can choose to be at x = 

D 2 [p] = -5, v V a {A{x, [p]))V a (A(x, \p})) + i [F<${x, [p]), } (55) 

where Fj£)(x, [p]) is related to the magnetic field by = e^xFj® and is given by 

#-sr/*"iST (56) 

Debye screening means that 

= ^ e- Mr /(Mr) (57) 

where the function f(x) has the property that for x » 1, f(x) ~ 1, and M is the mass gap 
related to the Debye length Id by M = y- 
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For a nonzero mass gap, the field outside the monopole core cannot be transformed to 
the radial form of (p6|). However close to the core and distances much smaller than the 
Debye length Id the field is close to the single monopole field outside the core and may be 
cast in the radial gauge. Furthermore at distances much larger than the Debye length, the 
field is close to zero and once again one may cast the gauge potential in the radial gauge 
trivially (i.e. with K(r) = 1). As mentioned in Section 4 the source of the even parity 
S-wave instability is the long range Coulomb field of the monopole. Since screening cuts off 
this Coulomb field and replaces it approximately by an exponential, one expects stability. 

The situation is in fact similar to that of the Yang-Mills-Higgs system in some respects. 
Recall that the S-wave stability of the 't Hooft-Polyakov monopole is ensured by the fact 
that the Higgs field rises exponentially to 1 beyond the core and cancels the negative \ tail 
of the gauge field, preventing the potential in the Schrodinger problem from being negative 
at large distances. In our problem the gauge field itself falls off exponentially to zero and 
thus the potential in ([36]) is positive at large distances. In this sense we have a dynamical 
Higgs effect. 

It is indeed possible to argue that the above argument for stability is self-consistent. Re- 
call that Debye screening means that the form of the magnetic field away from the monopole 



core is of the form |57|. We are unable to determine the precise form of the function f(x), 
but we can parametrize it in the following way. Inside the monopole core the magnetic field 
in the unitary gauge follows from the form of the single monopole vector potential 

B® = -£(1 - K (rA) 2 ) (58) 



The function K(x) is given by the solution in |16 |. Let us assume that this form of the 
magnetic field extends upto r = aX with some a > A. The magnetic potential far outside 
the core is an exponential, as required by Debye screening which gives 



B\ 6) = - — (1 + Mr) e~ Mr (59) 



where M is the mass gap generated in the plasma. We assume that this form of the magnetic 
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field extends from infinity upto the point r = aX. This gives the function K in the region 

r > aX 

K 2 {r) = 1 - (1 + Mr)e~ Mr r > aX (60) 

We further simplify the problem by replacing the function K(r) inside the core by an expo- 
nential so that we get 

K 2 {r) = e~ 2r/x r<a\ (61) 

The potential which appears in the s-wave stability operator is V(r) = 3K ^i~ 1 - This potential 
must be continuous at r = al which determines a implicitly through the equation 

e - 2a = 1 - (1 + MaX)e~ MaX (62) 

It may be now seen that the potential V(r) is negative in the region ^log 3 X < r < r 2 
where r 2 may be determined by the above considerations to be the solution of V(r) = for 
r > aX. This is approximately ^ 

We now test whether such a potential can have bound states. This may be done by 
applying the Bargmann criterion for absence of bound states 

JdrrV~{r)<l (63) 

where V~(r) stands for the negative part of the potential. A straightforward numerical 
integration then gives the result that the Bargmann bound is satisfied for a < 1.38. The 
fact that the critical value of a came out to be greater than one is an evidence for the self 
consistency of the picture. 

The above considerations are rough : we have made several simplifying assumptions. 
However these assumptions are dictated by the physics of the problem. The numbers quoted 
above are to be considered as indicative since they will change with different approximations 
to the function f(r) and K{r) inside the core. However the above calculation gives a self 
consistent argument in favor of the stability of monopoles in a neutral plasma. 
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G. Self consistent dynamical generation of monopole size 



As mentioned in Section the fugacity of the monopole gas depends on p and hence it has 
a dependence on the mass gap M and the cutoff A of the sine-Gordon theory. This means 
that the cutoff A is not independent, but determined in terms of M and g 2 , i.e. 

M = M(A, g 2 , z{\, M)) (64) 

implies that 

A = ~ 2 F(^) (65) 

g g 

where F is a function obtained by inverting (0). This equation now fixes the monopole core 
size A as a function of the mass gap M in a self-consistent manner. The above considerations 
must be considered as qualitative because the calculation of the function F in (|6~5D is beyond 
our present technology. 



V. DUAL REPRESENTATION AND THE DISORDER OPERATOR 

We now relate the euclidean formalism in terms of disorder operators introduced by 't 
Hooft Q . In the hamiltonian formalism the Schrodinger picture disorder operator &d(xo) is 
defined as an operator which creates a Z 2 magnetic vortex at the point xq in two dimensional 
space. More specifically it implements a singular gauge transformation which has the 
property that if we consider a closed spatial loop C parametrized by an angle 9 one has 

tt [xo] {6 + 2tt) = -Q [xo] {6) (66) 

when the loop C encloses the point xq. If C does not enclose x® the gauge transformation 
is single valued. Consider now the two point function of the Heisenberg picture disorder 
operators 

< & D {x)$ D (y) > (67) 
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Here x and y stand for three dimensional coordinates (including the euclidean time). This 
two point function is then a sum over all configfurations of the gauge fields which have a 
Dirac string singularity along a line joining x and y with a monopole of charge \ at the 
point y and an antimonopole of charge — | at the point x. It is crucial that the magnetic 
charges of these monopoles is half that of the monopoles which populate the vacuum. They 
have magnetic charges so that the Dirac string is visible by the lowest electrically charged 
quarks which couple to the gauge field. 

Repeating the steps which led to the sine-gordon representation with the difference that 
we have two external magnetic sources with charges ±| we easily get 

< & D {x)§ D {y) >=< e ^ x{y) ~ x{x)) > (68) 

the average on the right hand side in (|68D being performed in the sine-gordon theory. A 
similar identification holds for all higher point correlation functions of the disorder operators. 
Hence we can identify the disorder operator with 

$ D (a;) = e £x(«) (69) 

In fact the sine-gordon action may be now written in terms of $£> as 

S = 32^ / ^MlA* + M 2 ^d? + (& D ) 2 )} (70) 

upto an irrelevant constant. The field $d is not a conventional scalar field, since = 1. 
This non-linear Z 2 sigma model can be generalized to a linear sigma model by the addition 
of the term A / d 3 x(& D <& — l) 2 to the action (|70|). The action then exactly has the form 
conjectured in 0. 

The sine gordon theory thus is itself a dual representation of the original Yang-Mills 
theory. The action (f70"D has the global Z 2 symmetry $£) — > $£, which is spontaneously 
broken leading to magnetic disorder and confinement. This is simply the symmetry x — > — X 
of the sine-gordon model. It is clear from the action (|70|) that the dimensionless coupling 
constant is y -^. This is inversely related to the gauge coupling g as expected in a dual 
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formulation. The dual theory is weakly coupled when |f is small. In this limit the minima 
of the potential (cos x) break the Z2 symmetry spontaneously. 

Finally we note that the above construction of the disorder operators can be easily 
extended to SU (N) gauge theories following the treatment in |7| . 

VI. CONCLUSIONS 

We have argued that in 2+1 dimensional pure Yang-Mills theory Debye screening in a 
gas of regularized and dressed Wu-Yang monopoles provides a consistent picture of quark 
confinement. We have used the results of || that in a three dimensional Coulomb gas 
the charge density field always clusters, leading to Debye screening even for arbitrarily 
low temperatures. Our line of argument has been self-consistent in nature, because Debye 
screening in turn implies a screened magnetic field and hence the stability operator around 
a single dressed monopole is expected to have no negative eigenvalues. The mass gap thus 
obtained is non-perturbative and determines the monopolesize self-consistently. A related 
issue is that the mean configuration p[x) in the presence of a single monopole source is 
in general non-classical and hence the associated scalar potential x(x) does not satisfy a 
classical equation. Hence the explicit evaluation of the Wilson loops is not as easily done. 
However on general grounds the existence of a mass gap leads to qualitative conclusions 
that are similar to the case of the Yang-Mills-Higgs system. Finally we have obtained a 
representation of the disorder operators of the theory in terms of the sine-gordon field which 
leads to a dual representation of the gauge theory. 
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VII. APPENDIX I 



In this appendix we state the main results of |1 on the 3-dim. Coulomb gas. Theorem 
2.1 in Brydges ||, adapted to our notation states that: 

Given any c\ > 0, there exits C2 > such that for ^ < g 2 \ and z < \c\g % [z is the fugacity) 
the correlation functions of the density operator exits and clusters exponentially, i.e. there 
exist strictly positive constants M(z,g 2 , A), d = d(n') such that for m < n' 

< d exp (- 2 <„;%< n , \xi - Xj + a\ ■ M) 

M(z, g 2 , A) is the mass gap whose inverse is the Debye screening length. In the limit g 2 X — > 
oo one has the classical Debye-Huckel limit ^ — > 0. If we apply ( |7I|) to the 2-point function 
of the density, for separation of the order A, the monopole core size, which is the lattice 
spacing for the Coulomb gas we get 

|(Ap(0)Ap(A))|<cV MA (72) 

where Ap(x) = p(x) — {p(x)). ( ff2|) says that the fluctuations of p(x) are bounded and finite. 
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